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Specification 
Method of Designing Redox Flow Battery System 

5 Technical Field 

The present invention relates to a method of designing a redox flow 
battery system including a redox flow battery to force electrolytic solution 
to be fed to and discharged from its cells. Particularly, the present 
invention relates to a method of designing a redox flow battery system that 
10 can reduce a battery system loss. 

Background Art 

It is inherent in the electric power supply that power supply 
conforming to electricity consumption, what is called "simultaneous supply 
of equivalent electricity" is required. Meanwhile, developments of new 

15 energy resources, including power generation by wind and solar 

photovoltaic power generation, have been made increasingly in recent 
years. However, the power generation by wind, the solar photovoltaic 
power generation, and the like are all power sources which are irregular in 
output of power generation, due to which these power sources cannot serve 

20 as a satisfactory power supply source by themselves. Consideration is now 

being made on combination with some storage battery to stabilize the 
output of power generation. JP Laid-open (Unexamined) Patent 
Publication No. Sho 61*218070 describes the technique using a redox flow 
battery as the storage battery. 

25 This publication describes that a lead storage battery, which does not 
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need any moving parts, such as a pump and the like, is further annexed to 
the redox flow battery, in order to reduce a pump power loss and a loss 
caused by a shunt current, which contributes to reduction of combined 
efficiency of the redox flow battery, in response to input power. 

In general, the redox flow battery is used for equalization of load or 
for countermeasure to voltage sag. FIG. 8 shows an explanatory view 
showing an operating principle of a general redox flow secondary battery. 
This battery has a cell 100 which is separated into a positive electrode cell 
100A and a negative electrode cell 100B by a membrane 103 of an 
ion-exchange membrane. A positive electrode 104 and a negative 
electrode 105 are contained in the positive electrode cell 100A and the 
negative electrode cell 100B, respectively. A positive electrode tank 101 
for feeding and discharging positive electrolytic solution to and from the 
positive electrode cell 100A is connected to the positive electrode cell 100A 
through conduit pipes 106, 107. Similarly, a negative electrode tank 102 
for feeding and discharging negative electrolytic solution to and from the 
negative electrode cell 100B is connected to the negative electrode cell 100B 
through conduit pipes 109, 110. Aqueous solution containing ions that 
change in valence, such as vanadium ion, is used for the respective 
electrolytes and is circulated by using pumps 108, 111, to charge or 
discharge with an ionic valence change reaction at the positive and 
negative electrodes 104, 105. For example, when the electrolyte 
containing the vanadium ions is used, the following reactions occur in the 
cell during the charge or discharge of electricity- 

Positive electrode: V 4+ — » V 5+ + e- (Charge) V 4+ <— V 5+ + e- (Discharge) 



Negative electrode: V3+ + e - -> V 2+ (Charge) V* + + e - *- V 2+ (Discharge) 

However, the annex of the lead storage battery for the purpose of 
improving reduction of efficiency caused by weak power generation or load 
power at the time of electric charge or discharge or preventing increase of a 
system loss at the time of electric charge or discharge, as in the technique 
disclosed in the publication above, causes the problems of not only increase 
in production costs but also increase in scale of facilities. 

Driving the pump for feeding and discharging the electrolyte to and 
from the cells is absolutely necessary for the redox flow battery. Due to 
this, the conventional employs the annex of the lead storage battery to the 
power source which is irregular in output of power generation, while 
suffering from the disadvantages of increase in production costs and others. 

In general, a system loss means a total of a battery loss and a 
converter loss. The technique described in the publication cited above 
takes a pump power loss and a shunt current loss into consideration as 
factors of the battery loss, but takes no thought of an efficiency loss caused 
by a battery resistance (cell resistance). In view of this, the system that 
can provide further reduced loss is being desired. 

Meanwhile, there are some conventional facilities in which only the 
redox flow battery is combined with the power source which is irregular in 
output of power generation, without the lead storage battery annexed 
thereto, taking no thought of the reduction of efficiency caused by weak 
generated power or load power at the time of electric charge or discharge. 
In these facilities, no study is made on how great magnitude the storage 
battery should be for combination with the magnitude (output of power 



generation, variation in output of power generation, etc.) of the power 
source which is irregular in output of power generation, such as power 
generation by wind and solar photovoltaic power generation, and no design 
guideline thereof is provided. In the circumstances above, the redox flow 
battery is practically operated, combining with the storage battery of a 
reasonable magnitude, or specifically, a magnitude of output about one half 
the total output of power generation, but this operation induces a great 
system loss, leading to deterioration of system efficiency. 

Also, it is common that no study is made on how great the magnitude 
of the converter should be for combination with the magnitude (output of 
power generation, variation in output of power generation, etc.) of the 
power source which is irregular in output of power generation, such as 
power generation by wind and solar photovoltaic power generation, and no 
design guideline thereof is provided. For this, the redox flow battery 
combined with the power source irregular in output of power generation 
takes no thought of the converter loss. 

Further, no study has been made on the design guideline that can 
provide further reduction of the system loss for the redox flow battery used 
for equalization of load power or countermeasure to irregular variation in 
power consumption, such as voltage sag, as well as for the redox flow 
battery used for stabilization of output of power source which is irregular 
in output of power generation. 

Additionally, a method of designing an optimum redox flow battery 
considering not only the system loss but also reduction of production costs 
and scale of the facilities is being desired. 



It is a primary object of the present invention to provide a method of 
designing a redox flow battery that can provide a more optimal operation 
for the redox flow battery which is irregularly operated for stabilization of 
an output of power generation of a power source which is unstable in 
generated power, as well as for control of supply and demand of electricity. 

It is another object of the present invention to provide a method of 
designing a redox flow battery system that can reduce a system loss caused 
by weak generation power or load power at the time of electric charge or 
discharge, without using the lead storage battery. 
Disclosure of the Invention 

The present invention is defined on the basis of the knowledge given 
below. 

(1) Character to evaluate a system, such as a loss (loss characteristic), 
can be practically assumed as quadric functions. 

(2) External parameter pertinent to the character, such as a battery 
output for a smoothened output of power generation, has a probability 
distribution. 

(3) Expected value of the internal parameter pertinent to the 
character, such as a battery loss and a converter loss, is obtained from an 
average value of the internal parameter having the probability distribution, 
or standard deviation. 

(4) Optimal value of the internal parameter to optimize the character 
can be determined from the forgoing. To be specific, for example a 
specified output of battery or a specified output of converter for providing 
further reduction of the system loss is determined. The internal 



parameter for optimizing the character is expressed as n times of the 
standard deviation of the probability distribution. To be specific, for 
example the specified output of battery or the specified output of converter 
for providing further reduction of the system loss can be expressed as n 
times of the standard deviation of the probability distribution. 

The present invention is directed to a novel method of designing a 
redox flow battery system comprising a redox flow battery to force 
electrolytic solution to be fed to and discharged from its cells, comprising 
the steps of : first, determining an external parameter given by an 
operating condition of the redox flow battery a designer cannot choose 
voluntarily; then, determining an internal parameter given by a design 
condition of the redox flow battery the designer can design voluntarily; 
then, determining an average value of variables of the external parameter 
and standard deviation; and determining an optimum value of the internal 
parameter based on at least either of the resulting average value and the 
resulting standard deviation. 

It has been thought to be impossible hitherto to take an improvement 
in increase of loss caused by weak generation power or load power at the 
time of electric charge or discharge, without using the lead storage battery 
annexed to a power source which is irregular in output of power generation. 
Also, there has been no knowledge of an optimum method for designing a 
redox flow battery which is irregularly operated not only for stabilization of 
output of that power source but also for control of supply and demand of 
electricity. After having studied the possibilities of improvement from 
various angles, the inventors have found that the characteristics of the 



battery system above (reduction in system loss) can be optimized by 
determining magnitudes (specified output, number of cells, etc.) of the 
redox flow battery, magnitudes (specified output, etc.) of the DC/AC 
converter, such as an inverter, and other operating conditions (flow rate of 
electrolyte for each cell, temperature of electrolyte, etc.) by the specified 
method mentioned above, leading to the accomplishment of the present 
invention. In the following, the present invention will be explained in 
further detail. 

The external parameter is defined here as a parameter given by the 
operating condition of the redox flow battery the designer cannot choose 
voluntarily. The operating condition is properly varied depending on e.g. 
an atmosphere around the installation location of the redox flow battery, 
an intended application and purpose of the same battery, a user's demand, 
etc.. To be more specific, for example, an outside air temperature at the 
installation location of the redox flow battery is an example of the 
parameter the designer cannot choose voluntarily. The others include the 
smoothing of an output of power generation of the generating equipment, 
such as power generation by wind and solar photovoltaic power generation, 
which varies irregularly in output of power generation, the equalization of 
load, and the countermeasure to voltage sag. These are the parameters 
which are determined depending on the user's demand for the intended 
applications and purposes, and the power generation and the power 
consumption in those applications and purposes can hardly be determined 
by the designer. Accordingly, related parameters to the parameters cited 
above, such as, for example, the battery output for smoothing the power 



generation of the generating equipment that varies irregularly in output of 
power generation, and the battery output for smoothing the power 
consumption of the load that varies irregularly in power consumption, are 
also the parameters the designer cannot determine voluntarily. Further, 
kW capacity and kWh of the redox flow battery are the parameters which 
are determined by the user's demand and cannot be determined by the 
designer voluntarily. 

The external parameters that may be used include a parameter that 
takes a value that varies depending on a natural phenomenon, such as an 
outside air temperature at the installation location; a parameter that takes 
a value that depends on an irregular or stochastic phenomenon, such as, for 
example, a parameter that takes a value that varies temporarily, such as a 
battery output for smoothing irregular power generation or power 
consumption; and a parameter, such as kW capacity, that takes a fixed 
value. Thus, the external parameter may be only a parameter that takes 
a value that depends on a stochastic phenomenon, only a parameter that 
takes a fixed value, and, of course, a parameter that takes both of them, 
depending on the choice of the external parameter. Also, the external 
parameter may be used along or in combination of two or more. 

The present invention defines, in particular, the external parameters 
including at least one of the parameters that take values that depend on 
the stochastic phenomenon as recited above. The external parameter may 
additionally include a parameter that takes a fixed value. It is assumed 
that the parameter that takes a value that depends on the stochastic 
phenomenon is given an average value and standard deviation (or 
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variance). For example, when the battery output for smoothing the output 
of the power generation of the generating equipment that varies irregularly 
in output of power generation or the battery output for smoothing the 
power consumption of the load that varies irregularly in power 
consumption is chosen as the external parameter, the average value and 
the standard deviation can be determined from individual outputs that 
vary temporarily. The term of "individual outputs that vary temporarily" 
is intended to include, for example, a battery output (kW) at a point in time 
in a date, a sequent battery output (kW) at + one second after that, a 
further sequent battery output (kW) another + one second after, ... . For 
example, when an outside air temperature at the installation location is 
chosen as the external parameter, the average and the standard deviation 
can be determined from individual ambient temperatures that vary 
temporarily. The term of "individual ambient temperatures that vary 
temporarily" is intended to include, for example, a temperature (K) at a 
point in time in a date, a sequent temperature (K) at + one second after 
that, a further sequent temperature (K) another + one second after, ... , 
when examined at K (Kelvin) by the minute. 

The external parameter having the average and the standard 
deviation as chosen for the application and purpose include, for example, 
those cited below. A battery output for smoothing the output of power 
generation of the generating equipment which varies irregularly in output 
of power generation can be cited as an external parameter for the 
application and purpose of smoothing the output of power generation of the 
generating equipment, such as power generation by wind and solar 



photovoltaic power generation, which varies irregularly in output of power 
generation. A battery output for smoothing the power consumption of the 
load that varies irregularly in power consumption can be cited as an 
external parameter for the application and purpose of equalizing the load 
power. An outside air temperature at the installation location can be cited 
as an external parameter for the application and purpose of 
countermeasure to voltage sag. 

The phase of "smoothing the output of power generation" used herein 
is intended to mean that when an output of power generation exceeds a 
threshold as is preset for the output of power generation, the surplus 
output exceeding the threshold is charged in the battery, while on the other 
hand, when an output of power generation is less than the threshold, the 
output corresponding to the shortage is discharged from the battery. Also, 
the phase of "smoothing the power consumption" used herein is intended to 
mean that when power consumption exceeds a threshold as is also preset 
for the power consumption, the output corresponding to the shortage 
caused by the power consumption exceeding the threshold is discharged 
from the battery, while on the other hand, when power consumption is less 
than the threshold, the surplus output is charged in the battery. The 
same or different thresholds may be used for charging and discharging the 
battery. Also, the thresholds may be varied depending on the time 
required for the power generation and the output situation thereof. 

The internal parameter is defined here as a parameter given by the 
design condition of the redox flow battery the designer can design 
voluntarily. To be more specific, the internal parameters that may be 

10 



used include, for example, a specified output of the battery, the number of 
batteries, the number of cells, a fluid volume of an electrolytic solution 
reservoir tank, a flow rate of the electrolytic solution for each cell, a 
temperature of the electrolytic solution, a specified output of a DC/AC 
5 converter for converting the battery output, and the number of DC/AC 

converters for converting the battery output. Thus, the internal 
parameter is a parameter that takes a fixed value, differently from the 
external parameter. The internal parameter may be used alone or in 
combination of two or more. 
10 Now, a method of calculating an optimum value of the internal 

parameter is explained concretely. 

(1) The case where the external parameter is only a parameter that takes 
a fixed value- 

Let the external parameter be represented as xk X2n • xi N X2^ "* are 
15 all taken as fixed values different from each other- for example, xi : kW 

capacity, X2 : kWh capacity, 

Let the internal parameter be represented as yi N y2^ •'*. y^ y2> *•• are 
taken as the parameters different from each other- for example, yi : a 
specified output of battery, y2 : a specified output of a DC/AV converter, 
20 Let the characteristic function to evaluate the system be represented 

as f(xu X2x yi> y2^ •••). f includes, for example, a system loss (efficiency), 
cost, and size, f includes at least one of these parameters. 

Then, f is optimized when f takes a local extreme value, i.e., when d fl 
dyi = 0 or when yi takes a maximum value or a minimum value in a 
25 variable range. In the latter case, a value at a boundary of the variable 
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range can be cited as a possible maximum value or a possible minimum 
value. 

Hence, when the external parameter is a parameter that takes a fixed 
value, the optimum value of f can be easily determined by partial 
differentiation as mentioned above. 

<2> The case where the external parameter includes a parameter that 
takes a value that depends on an irregular or stochastic phenomenon- 
Let the external parameter be represented as xk X2^ * *Xi N It is 
assumed that xi N X2> * * xu ••• are parameters different from each other. 
Some of these parameters Xi are the parameters that take values that 
depend on a stochastic phenomenon. For example, xki(i = kl) is the 
battery output for smoothing output of power generation of the generating 
equipment that varies irregularly in output of power generation, and Xk2(i 
= k2) is the outside air temperature at the installation location of the 
battery, etc.. Also, Xi is the parameter to which the average and the 
standard deviation are given. The average is represented as Xi ave» and the 
standard deviation is represented as o X i(variance o X i 2 ). The other external 
parameters xk X2n *** are the parameters which take fixed values. 

Let the internal parameter be represented as yu y2x ** in the same 
manner as in the above. 

Then, if the characteristic function to evaluate the system f(xi N X2> 
yu y2> "•) can be rewritten in the form of a new characteristic function g(xu 
X2> • > xkiave. a xki 2 ^ xk2aveN o xk2 2 , " ^ *iave, o x i 2 , yu y2^ " ) considering 
the average Xi ave and the variance o xi of the external parameter Xi, then g 
can be optimized when g takes a local extreme value, when d f/ d yi = 0, or 
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when yi takes a maximum value or a minimum value in a variable range. 

Consequently, a possible condition for allowing f(xi N X2s • ^ y;u y2^ 
"0 to be rewritten in the form g(xi N X2s • > Xkiave> o xki 2 ^ xk2ave> Oxk2 2 > 
• Xiave^ a x i 2 > yu y2^ ' ' 0 is set up, first. 

Here, the parameter xiof the external parameters is assumed to have 
a value that depends on an irregular or stochastic phenomenon and also 
have an average value Xiave and a variance o xi • Then, it is assumed that 
Xi can obtain individual temporarily variable values Xii, Xi2, xin. 

Then, it is assumed that the average Xi ave can be given by Eq. 1.1 and 
the variance o X i 2 can be given by Eq. 1.2. 

x +s + .■■+*,, 
x, ove = ^ • • Eq- 1- 1 

i O, a™ - f + (X ove - x, 2 ) 2 + . . . + (x, ave - x, ) 2 

O v = : ... rLQA.Z 

N 

When the individual temporarily variable values Xii N Xi2x ■ ^ XiN are 
rewritten in the form of Xii=xiave+ Ak Xi2 — Xi ave+ A 2^ XiN— Xiave+ An 
using an average Xi ave and the difference A j from the average Xi ave? they 
can be rewritten as the following Eq. 1.3 and Eq. 1.4. 

From Eq. 1.1, A , + A 2 + ■— + A N =0 ...Eq. 1.3 



(A,) 2 +(A 2 y+... + (A N ) 2 
N 



From Eq. 1.2, K ~" ' "~ 2 \' "' ' '~ N/ = a x 2 ...Eq.1.4 



Then, when the characteristic function f(xi s X2^ yi> ■••) to 
evaluate the system is expanded around Xiave, or with respect to Xi ave + A, 
using Taylor expansion, Eq. 1.5 is given. 
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f(x,,x 2 ,...,x i aw +A,y,,y 2 ,...) = 

f(x„X 2 ,..., Xi „, yi ,y,,..^ — .4 +1 ^.j + -—.j + ... 

...Eq.1.5 

When expected values f a ve of the characteristic function f with respect 
to the individual values xik Xi2^ xin of the external parameter Xi are 
determined by using Eq. 1.5 above, Eq. 1.6 is obtained. 
1 ^ 



fa™ =-jj Y t f{x tt x 2t ... t x ljt y lt y 2t ...) 

= f( X \ i X 2>' — > X i ave > J^l ' > • * •) 



df ( A ]+ A 2 +... + A N \ 1 3 2 f 

dxA N ) 2dx, 



(A,) 2 +(A 2 ) 2 +... + (*„) 



2 \ 



+ 



1 d 3 f / (A,) 3 +(A 2 ) 3 +... + (A„) 



6 dx 



2 

i V 
3 \ 



N 



N 



... EqA.6 

When Eq. 1.3 and Eq. 1.4 are substituted for Eq. 1.6 above, Eq. 1.7 is 
obtained. 

fave ~ f( X ] 9 X 29'"9 X , ave > .Vl 5 ^2 ' * * •) 

, 1 tL. a 2 , i a 3 / r (A,) 3 +(A 2 ) 3 +... + (A A ,) 3 > | | ^ 



2 dx i 



6dx, 



N 



... EqA.l 

Hence, it is derived from Eq. 1.7 that when third or more order 
derivatives of f(xu X2s •xu yu y2> • ) at xi becomes zero, i.e., when the 
following Eq. 1.8 is satisfied, it is made possible to rewrite f(xi> X2> •> Xi > 



yu y2^ • ■) in the form of g(xu X2^ • Xiaves o xi 2 . • ^ yu y2. 



•). 



dc/ dc; 

This is the case where f(xi N X2^ xk yu y2> •••) is in the form of a 
quadratic equation with respect to Xi. Specifically, when 
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f(xu X2> • • • % xi, yi. ...) 

=a(xu x 2 , yis y 2 ^ ...) *Xj 2 

+b(xj, x 2 , yi. y 2 , ...) -Xj 

+c(xi, x 2 , .... yi, y2s -J (Eg. 1.9), 

g(xi X2^ .... Xiave^ O x i 2 s .... yi, y2 > ...) 

= a-Xi 2 +b -Xi +c +a -a J (Eq. J JO). 

This indicates that the condition for allowing f(xi N X2. yu y2. ") 
to be rewritten in the form of g(xu X2. * •> Xkiave> o x ki 2 > Xk2aveN o x k2 2 > 
Xiave> axi 2 . yi. y2. * ") is obtained when f(xu X2. * ^ yi, y2. ••*) comes 
to be in the form of quadratic equation with the external parameter xi 
having the average and the variance. 

It should be noted that if f(xu X2> xk yi. y2^ • ) does not exactly 
comes to be in the form of quadratic equation with respect to Xi but comes 
near the quadratic equation, then second order derivatives will also be 
taken into consideration. Practically, there seems to be many cases where 
the quadratic equation is sufficient to approximate the actual properties. 
But, it is derived from Eq. 1.7 that a cubic or more equation can also be 
used to approximate the actual properties, considering the effect of terms 
higher than a quadratic term, specifically : 

(A 1 ) 3 + (A 2 ) 3 +... + (A A/ ) 3 (A ] ) 4 +(A 2 ) 4 +... + (A N ) 4 
TV ' N 

To give actual examples in the following- 

(l) The case where f comes to be in the form of linear expression with 
respect to xi- 

f=a-Xi + b — » g=a-xiave + b 
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Hence, g has an optimum value obtained by the average Xiave. 

(2) The case where f comes to be in the form of quadratic expression with 
respect to Xi- 

f=a-Xi 2 + b*Xi + c — > g— a-Xi ave 2 + b -Xi ave + c+ a- a xi 2 
Hence, g has an optimum value obtained by the average Xi ave and the 
variance a X i 2 . 

(3) The case where f comes to be in the form of cubic expression with 
respect to x^ 

f = a • Xi 3 + b * Xi 2 + c • xi + d — > 



f = a m x l ave +3a • x,. ove - a xi +a • — ' 2 — = — ■ 

+ b ' x, a J +b - a J +cx, aw + d 

Hence, g has an optimum value obtained by the average Xi ave, 



the variance o x i 2 , and ■ 



N 

Although the case where a single external parameter having a 
value that depends on an irregular or stochastic phenomenon or a 
single parameter to which average and standard deviation are given is 
chosen among the external parameters has been described above, there 
may be the case where two or more external parameters are chosen 
among the external parameters. In this case, f can be rewritten in 
the form of quadratic expression concerning those parameters. Take 
the case where two external parameters to which average and 
standard deviation are given are chosen, for instance. 

Let the external parameters to which average and standard 
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deviation are given be represented as x, y. Also, let their averages 
and variants be represented as Xave> Vave\ o x 2 n o y 2 , respectively. It is 
assumed that these parameters x, y are not correlated with each other. 

Then, assuming that the characteristic function f(x, y) can be 
arranged in the form of quadratic function with two unknowns with respect 
to x and y, and letting a coefficient be represented as aij, it can be rewritten 
as follows. 

f=a22 # x 2 -y 2 + a2i # x 2 -y 1 + a2o-x 2 -y° 
+ ai2 • x 1 • y 2 + ai i • x 1 • y 1 + aio • x 1 • y° 
+ ao2 • x° * y 2 + aoi ■ x° • y 1 + aoo • x° * y° 

When rearranging this in the form of x = x a ve+ A X n y — yave+ A y by 
using difference A between the individual temporarily variable values 
and the average and then substituting them for the function f, 

f = a22 • (Xave + A x ) 2 ' (yave + A y ) 2 + a21 ' (x av e + A x ) 2 * (yave + A y) 
+ a20-(Xave+ A x ) 2 

+ ai2 * (Xave + A x) * (yave + A y ) 2 + an • (x a ve + A x ) • (y a ve + A y ) 
+ aiO-(x aV e+ Ax) 

+ ao2 • (yave + A y ) 2 + aoi * (y ave + A y ) + aoo 

When the terms of the function f above are expanded by using Taylor 
expansion, consideration of the effects of the averages and variants of the 
external parameters must be given to the expanded terms including A x ^ A 
yN A A A y 2 . A x A y> A x Ay 2 , A x 2 A y . A x 2 A y 2 as the coefficients. 

The following change of variables is given to the respective expanded 
terms of the function f. 

® Ax-^O Ay^O Ax 2 ^<Tx 2 Ay 2 ^ay 2 AxAy^O 
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© A x Ay 2 ^0 Ax 2 A y ^0 Ax 2 A y 2 -*a x 2 a y 2 

In the change of variables CD , A x ^(h A y — >0, A x 2 — > o A A y 2 — > o y 2 
are determined from the definitions of the averages and the variances. A x 
A y ^0 is determined from the definition of "non-correlated". As to the 
5 change of variables (2), it is assumed that the same change as above is 

given. 

f(x, y) may be rewritten in the form of g(x ave ^ yave^ o x 2 > o y 2 ) by using 
the change of variables of ® and (2) above. 

Following the sequences above, the calculation is made for 

10 optimization of the system. Particularly, when at least one of a specified 

output of battery, number of batteries, a specified output of DC/AC 
converter for converting the battery output, and number of DC/AC 
converters for converting the battery output is obtained from the battery 
output of the redox flow battery for smoothing output of power generation 

15 of generating equipment which varies irregularly in output of power 

generation, the following procedures may be taken. First, standard 
deviation of output distribution of the redox flow battery used for 
smoothing output of power generation of the generating equipment which 
varies irregularly in output of power generation is determined. Then, 

20 based on the resulting standard deviation, the at least one of the specified 

output of battery, the number of batteries, the specified output of DC/AC 
converter, and the number of DC/AC converters is determined. 

It is preferable that the specified output of the DC/AC converter is in 
the range of not less than 1 time to not more than 4 times, or preferably in 

25 the range of not less than 1 time to not more than 2.5 times, of the 
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standard deviation of the output distribution of the battery with respect to 
the smoothed output of power generation. Also, it is preferable that the 
specified output of the battery is in the range of not less than 0.7 time to 
not more than 2 times, or preferably in the range of not, less than 0.9 time 
5 to not more than 1.5 times, of the standard deviation above. When the 

specified output of the converter or that of the redox flow battery satisfies 
the value specified above, the system loss can be reduced further, leading to 
further improvement in system efficiency. The specified output of the 
battery can be changed by changing area or dimension of the electrode, 

10 changing the number of cells arranged in series/parallel, and so on. The 

specified output of the converter can be changed by changing capacity of a 
semiconductor device used in the converter, changing the number of 
devices arranged in series/parallel, and so on. 
Brief Description of the Drawings 

15 FIG. 1 is a graph showing a relation between a loss characteristic of a 

redox flow battery and a loss rate of the same; 

FIG. 2 is a histogram of an output of a battery of a battery system; 
FIG. 3 is a graph showing a loss characteristic of a converter with 
respect to an output of an AC end-cell battery; 

20 FIG. 4 is a graph showing a relation between a specified output of the 

converter and a loss characteristic of the same; 

FIG. 5 is a graph showing the loss characteristic of the redox flow 
battery with respect to the output of the AC end-cell battery; 

FIG. 6 is a graph showing a relation between a coefficient of an 

25 approximation formula and a flow rate of electrolyte per unit time for each 
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cell; 

FIG. 7 is a graph showing a relation between a specified output of the 
redox flow battery and the loss characteristic of the same; 

FIG. 8 is a schematic view showing an operating principle of the redox 
5 flow battery. 

Best Mode for Carrying out the Invention 

In the following, an embodiment of the present invention is described. 
(Outline of a calculation method of an optimum value of internal 
parameter) 

10 (l) Choose xi, X2 as external parameter. Choose yi, y2 as internal 

parameter. 

It is assumed that the external parameter xi depends on an irregular 
or stochastic phenomenon and an average value xiave and a variance a x i 2 
are given. It is also assumed that the external parameter X2 is given a 
15 fixed value. 

(2) A characteristic function to evaluate the system f (xi, X2, yi, y2) is 
determined and f is arranged in the form of quadratic equation with 
respect to xi. 

f=a(x2^ yu y2)xi 2 + b(x2 N yi. y2)xi + c(x2> yu y2) 
20 (3) A new characteristic function g considering an average value xiave 

and a variance a x i 2 of the external parameter xi is given as follows. 

g=a(x2> yu y2)xiave 2 + b(x2. yu y2)xi ave + c(x2. yu y2> + a(x 2 > yu 
y2) a xi 2 

(4) yiwhen the new characteristic function g takes an extreme value, 
25 i.e., yi when d f/ d yi = 0, or a maximum value or a minimum value in a 
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variable range of yi is a candidate for optimizing g. Then, yi to optimize g 
is an optimum value of the internal parameter. 

In the following, the external parameter, the internal parameter, and 
the characteristic function f are respectively explained concretely. In the 
5 example, reference is made to a calculation method of an optimum value of 
the internal parameter that can optimize the system efficiency for the 
purpose of smoothing of an output of power generation of the generating 
equipment, such as the power generation by wind and the solar 
photovoltaic power generation, which varies irregularly in output of power 
10 generation. 

(The case where an output of battery is constant to a smoothed output of 
electricity) 

An output of battery (hereinafter it is referred to as battery output) 
for the smoothed output of power generation is chosen as the external 

15 parameter. Here, the battery output is considered to be constant. Also, a 

specified output of the redox flow battery, a specified output of the 
converter, a flow rate of electrolyte, and a temperature of the electrolyte 
are taken as the internal parameter. 

With the system efficiency at the time of electric charge or discharge 

20 defined by Q battery output/(battery output + loss) and (2) (battery output 

— loss)/battery output, and a loss rate defined by (3) loss/battery output, 
the condition for maximizing the system efficiency or the condition for 
minimizing the loss rate is determined with respect to each definition. In 
any case, the system efficiency and the loss rate are partial-differentiated 

25 with respect to the battery output and the internal parameter to determine 
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the condition for maximizing the system efficiency or the condition for 
minimizing the loss rate. ® is described below, first. 

Let the battery output be represented as x, the internal parameter as 
y, the loss as fix, y), and the system efficiency at the time of electricity 
charge or discharge as tj = x/(x + f). When the system efficiency is 
maximized, the efficiency takes a local extreme value. Hence, the 
condition for maximizing the system efficiency is d 77 / d x = 0, d r\ I d y = 0, 
or leads to the following equations 2.1 and 2.2- 

~ l-(x + /)-x-(l + -^) 

a ~ / , ~ u Eq. 2.1 

dx (x + /) * 



0. (*+/)-*. (0 + 7 A 

^= ^-=0 • • • Eg. 2.2 

dy (x-f/) 2 



It follows from the denominator (x + f) 2 >0 in Eq. 2.1 and Eq. 2.2 that 
the numerator = 0. Hence, with respect to (J) , the condition for 
maximizing the system efficiency is given when satisfying the following 
equations 2.3 and 2.4. 

^— = 0 • • • Eg. 2.3 

x dx 

^=0 ...Eq. 2A 

dy 

The same calculation was made to determine the conditions with 
respect to (2) and (3). It was found therefrom that the conditions derived 
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from the fact that when the system efficiency is maximized, the efficiency 
takes a local extreme value, or when the loss rate is minimized, the loss 
rate takes a local extreme value are both the same as the condition above. 
To be specific, those conditions also satisfied Eq. 2.3 and Eq. 2.4. This 
means that even when the definitions of the system efficiency are different 
and the figures themselves of the efficiency are different, or even when 
considering from the viewpoint of the loss rate, the system efficiency is 
maximized (the loss rate is minimized) with the same battery output and 
internal parameter. Hence, the loss rate can be determined by doing a 
mathematical simulation of the relation between the loss and the output 
and internal parameter of the redox flow battery and then dividing the 
resulting loss by the output of the AC end-cell battery. Shown in FIG. 1 is 
an example of a graph showing a relation between a loss characteristic of 
the redox flow battery and a loss rate of the same. 

As shown in FIG. 1, when electricity is charged (when output of the 
system is positive), the loss rate of the battery and that of the converter are 
minimized when the output of the system is little under 150%, from which 
it can be seen that the battery loss and the converter loss can be minimized 
at that time. When electricity is discharged, the loss rate of the battery 
and that of the converter are minimized when the output of the system is 
about ~100%, from which it can be seen that the battery loss and the 
converter loss can be minimized at that time. 

From the foregoing, it can be seen that in the case where the battery 
output is constant, the battery output and the internal parameter that can 
allow maximization of the system efficiency (minimization of the loss rate) 
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or can allow minimization of the system loss can be determined. 

(The case where the battery output depends on a stochastic phenomenon) 

Next, consideration will be given to the case where the output of 
power generation varies irregularly, or specifically, the case where the 
battery output has a probability distribution. Let the battery output at a 
regular time interval (At) be represented as xu X2> ***> x n , the internal 
parameter as y, the loss as f(xi,y),f(x2,y),...,f(xn,y), and the system efficiency 
at the time of charge or discharge of electricity as 77 . The battery output 
xi,X2,..., x n has a probability distribution, so the average x ave gives Eq. 3.1 
and the variance a 2 gives Eq. 3.2. Then, when f(xi,y) is expanded around 
Xave, using Taylor expansion, Eq. 3.3 is given. 

x aw =^f- ...Eg. 3.1 



(T 2 = E(X, X ove ) 2 Eg32 

n 



/(x^y) = f{x avc y) + — (x, - x ahc ) + — - * flVC ) + Eg.33 

ox lux 

It follows from Eq. 3.1 that Zxi—n 'x ave (Eq. 3.4) and from Eq. 3.3 that 
EfCxi.y) leads to Eq. 3.5. 

T.f{x u y)= n.(f(x avc y) + -?^ T a 2 +...) • • ■ Eq .3.5 

2 ox 

Here, if n * Xave of Eq. 3.4 is taken as x in the case where the battery 
output is constant and the corresponding part of Eq. 3.5 to Eq. 3.6 given 
below is taken as f(x,y), then the case where the battery output has a 
probability distribution can also be dealt with in the same manner as in the 
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case where the battery output is constant. 

»-(/(W) + ^4-°- 2 +-) -Eq.3.6 
Lux 

Then, if f(x,y) is taken as a function with x, and the loss characteristic 
of the battery and that of the converter are practically taken as quadratic 
5 functions, then third or more order derivatives can be taken as 

substantially zero. To be specific, in the case where the probability 
distribution is considered, the loss characteristic which is to be dealt with 
by f(x)=a-x 2 +b'x+c (Eq. 3. Z^with respect to a constant battery output x 
can be dealt with by g(x a ve, a 2 )=a *x 3ve 2 +b -x ave +c +a * a 2 (Eq. 3.8) with 
10 respect to the battery output's average x ave and variance o 2 . 

(Example) 

A redox flow battery system including generating equipment that 
varies irregularly in output of power generation, a redox flow battery for 
smoothing the irregular output of power generation, and a DC/AC 

15 converter for converting the battery output was produced, and the loss 

characteristic of the battery and that of the converter were examined. 

A power generation by wind having a maximum output of power 
generation of 400kW (specified output of power generation' 275kW) was 
used as the generating equipment. A redox flow battery having the 

20 construction shown in FIG. 8 was produced as the redox flow battery and 

set at a specified battery output of 170kW (about 60% of the specified 
output of power generation). A DC/AC converter having a specified 
converter output of 275kW was used. A histogram of 8-hour battery 
output for the battery system comprising the generating equipment, the 
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redox flow battery and the converter is shown in FIG. 2. The distribution 
characteristic is shown below. The standard deviation was determined by 
smoothing the output of power generation of the generating equipment by 
using the battery, followed by the output distribution of the battery with 
respect to the resulting smoothed output of power generation. The term of 
"specified output of battery" means an output at which the system 
efficiency is maximized during the load smoothing operation of the battery. 
Also, the term of "specified output of converter" means an output at which 
the conversion efficiency of DC— >AC or AC— ^DC is maximized. 

Distribution characteristic of charged/discharged power at the time of 
smoothing of output of power generation* 

Smoothed output of power generation by wind* 200kW 
Average value of battery output • 0.9kW 

Standard derivation of battery output * 54.1kW 

Average of absolute value of battery output : 44.4kW 
(Determination of number of converters) 

Then, the loss characteristic of the DC/AC converter was determined 
(by factory test). The results are shown in FIG. 3. It is found from an 
approximated curve shown in FIG. 3 that it follows from Eq. 3.7 that the 
converter loss f(x)= 1.22 X 10" 4 -x 2 - 1.74X 10*g- x + 4.58 with respect to the 
constant battery output x. Hence, it follows from Eq. 3.8 that an expected 
value of the loss characteristic with respect to the battery output at the 
time of the smoothing operation of the power generation by wind having 
the probability distribution g(x ave , o 2 ) = 1.22 X 1CM • (0.9 2 + 54. 1 2 ) - 1.74 X 
10-6-0.9 + 4.58 = 4.94kW (x a ve = 0.9, a 2 = 54.1 2 ). It is to be noted that in 
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g(xave, or 2 ) of the equation above, the average value x ave has negligible 
effects on the loss. 

The internal parameter that can keep the loss characteristic of the 
converter of the specified output of 275kW is discussed further concretely. 
For example, the case where n converters of the specified output of 275kW 
are arranged in parallel is considered as the internal parameter. Then, 
the output for each converter is 1/n. Then, the expected value of the loss 
characteristic for each converter is f(x) = 1.22 X 10" 4 -(x/n) 2 — 1.74 X 10<Kx/n) 
+ 4.58. Then, the expected value of the total loss characteristics in the 
case where n converters of the specified output of 275kW are arranged in 
parallel is f(x,n) = 1.22X 10" 4 /n-x 2 -1.74X 10"6-x + 4.58-n. Hence, g(x ave , a 
2 ,n) = 1.22 X 10" 4 /n • (x ave 2 + a 2 )- 1.74 X 10 6 • x ave + 4.58 • n. It follows from 
average: x ave = 0.9, variance a 2 = 54.1 2 that g(n) = 1.22X 10-4/n-(0.9 2 + 54.1 2 ) 
- 1.74 X 10" 6 • 0.9 + 4.58 • n = 0.3572/n - 1.57 X 10 6 + 4.58 • n. Then, when 
finding n such that d g/3n-0, -0.3572/n 2 + 4.58 = 0, so n = 0.28. 

A variable range of n is 0^n^°°. Then, 0 can be cited as a possible 
minimum value n, and infinity can be cited as a possible maximum value, 
but when n— >0, g^°°, and when n^°°, g— so n^0.28 cited above is an 
optimum value. 

This means that when 0.28 converter of the specified output of 275kW 
is used, the loss can be minimized. In other words, when a converter of a 
specified output of 275 X 0.28^ 77kW is used, the loss can be minimized. 
FIG. 4 shows the expected value of the loss characteristic of the converter 
(converter loss). 

It can be seen from FIG. 4 that in this example, when the specified 
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output of the DC/AC converter is in the range of about 50kW to about 
200kW (about 1 to about 4 times of the standard deviation of 54.1kW), or 
particularly in the range of about 50k W to about 120kW (about 1 to about 
2.2 times of the same), the converter loss is reduced and thus the system 
loss can be reduced. From the foregoing, it can be seen that when the 
specified output is set at 77kW, the system loss can be minimized. 
(Determination of number of batteries and flow rate per unit time) 

Then, the loss characteristic of the redox flow battery was determined. 
The results are shown in FIG. 5. Since the loss characteristic of the 
battery varies depending on a flow rate of electrolyte per unit time for each 
cell, the loss characteristic of the battery (battery loss) is determined, 
varying the flow rate as follows. 

® 0.80 liter/min. -cell 

® 1.00 liter/min. -cell 

(3) 1.15 liter/min. -cell 

The loss characteristic of the battery can be determined from an 
approximated curve shown in FIG. 5 in the same manner as in that of the 
converter mentioned above. Take (2) the flow rate of 1.00 liter/min. • cell 
for instance, it follows from Eq. 3.8 that g(x ave , a 2 ) = 4.33 X 10' 4 -(x a ve 2 + a 2 ) 
— 2.28 X 10' 2 -Xave+ 11.6 with respect to the battery output at the time of 
smoothing operation of the power generation by wind having the 
probability distribution. Hence, it follows from x aV e — 0.9, a 2 = 54.1 2 that 
the expected value of the loss characteristic g(x ave , a 2 ) = 4.33X 10" 4 -(0.9 2 + 
54.1 2 )-2.28X10- 2 -0.9+11.6 = 12.8kW. Similarly, in the case of® 0.80 
liter/min. -cell, it follows from Eq. 3.8 that g(x ave , cr 2 ) = 4.72 X 10' 4 -(x ave 2 + o 
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2 ) — 2.35 X 10" 2 -Xave + 9.2, so the expected value of the loss characteristic 
g(xave,a 2 ) = 10.6kW. 

The internal parameter that can keep the loss characteristic of the 
battery of the specified output of 170kW is discussed further concretely. 
For example, the case where n batteries of the specified output of 170kW 
are arranged in parallel is considered as the internal parameter. Also, the 
flow rate of electrolyte per unit time for each cell is also determined for the 
operating condition. Then, the output for each battery is 1/n. Then, the 
loss characteristic for each battery is f(x) = a(x/n) 2 + b(x/n) + c. Then, the 
expected value of the total loss characteristics in the case where n batteries 
are arranged in parallel is f(x,n) = ax 2 /n + bx + cn. Coefficients a, b, c of 
f(x,n) are derived from an approximate expression shown in FIG. 5. Also, 
the relation between the coefficients a, b, c and the flow rate of electrolyte 
per unit time for each cell is shown in FIG. 6. 

The graph of FIG. 6 showing the coefficients is plotted, using the 
coefficients of their respective quadratic functions shown in FIG. 5. Take 
the coefficient a, for example- when the approximation formula is 
determined from the quadratic coefficients of 4.14, 4.33, and 4.72 of their 
respective quadratic functions CD, (2) and (3) shown in FIG. 5, it can be 
determined to be a( X 10' 4 ) = (-1.673) • a' + 6.0418, using the variance a'. 
The same applies to the coefficients b and c. 

Then, it follows from FIGS. 5 and 6 that f(x,n,L) = {(-1.673L + 6.0418) 
X 10- 4 }/n • x 2 + {(-0.073L + 2.3918) X 10" 2 } • x + {(1.3622L - 0.1795) X 10 1 } • n, 
where L(liter/min.) is a flow rate of the electrolyte for each cell. 

Hence, g(x av e, o 2 ,n,L) = {(-1.673L + 6.0418) X 10" 4 }/n ■ (x ave 2 + a 2 ) + 
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{(-0.073L + 2.3918) x 10 2 } • x aV e + {(1.3622L - 0.1795) X 10 1 } ■ n. It follows 
from average x ave = 0.9, variance a 2 = 54. 1 2 that g(n,L) = (-0.490L+ 1.769)/n 
+ (-0.0657L + 2.153) X 10 2 + {(1.3622L - 0.1795) X 10 1 } • n. Then, when 
finding n such that 3g/<9n=0, d g/ d n=-(-0.490L+1.769)/n 2 +(l3.622L- 
5 1.795) = 0, so, n = ^{(-0.490L+1.769)/(l3.622L- 1.795)}. 

Also, when finding n such that d g/ 3 L = 0, d g/ d L=-0.490/n — 
0.000657 + 13.622-n = 0. Then, it follows from n>0 that n=0.19. Then, 
it follows from 0.19 = J~ {(-0.490L + 1.769)/(13.622L-1.795)} that L = 
1.87(liter/min.). 

10 Here, if the flow rate of electrolyte of the redox flow battery shown in 

this example is assumed to have a possible range of 0.8 to 1.15 liter/min., 
then L^1.87 is a value outside the possible range. Hence, a possible 
optimum vale in the variable range of L is- 

Combination of n such that d g/ d n = 0 and the minimum value of 0.8 

15 in the variable range of L, i.e. (n, L) = (0.39, 0.8), or 

Combination of n such that <9g/3n=0 and the maximum value of 
1.15 in the variable range of L, i.e. (n, L) = (0.29, 1.15). 

By substituting the value above for g(x a ve,o 2 ,n,L), an optimum value is 
obtained when (n,L) = (0.39,0.8). 

20 This means that when 0.39 battery of the specified output of 170kW is 

used and the flow rate of electrolyte is set at 0.8 liter/min., the loss can be 
minimized. In other words, when a battery of a specified output of 170 X 
0. 39 ^66. 3k W is used, the loss can be minimized. FIG. 7 shows a graph of 
expected values g(x ave , a 2 ,n,L) (by letting x ave = 0.9 and valiance o 2 = 54.1 2 ) 

25 when the flow rate of the electrolyte is taken to be in a possible range of 0.8 
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to 1.15 liter/min.. 

It can be seen from FIG. 7 that in this example, when the flow rate 
shown in (2) above is 1.00 liter/min. • cell, the specified battery output is in 
the range of about 40kW to about lOOkW (about 0.74 to about 2 times of 
the standard deviation of 54.1kW), or particularly in the range of about 
50kW to about 80kW (about 0.92 to about 1.5 times of the same), the 
battery loss is reduced and thus the system loss can be reduced. From the 
foregoing, it can be seen that when the flow rate is set at 0.8 liter/min. • cell 
and the specified battery output is set at 66.3kW, the system loss can be 
minimized. 

From the foregoing, it was found that when the system charges or 
discharges electricity irregularly, with the average value of the 
charged/discharged power of the system of substantially zero, for 
stabilization of the power generation by wind and the like, the design of the 
internal parameters, such as the specified output of the battery, the 
specified output of the DC/AC converter, and the flow rate of the electrolyte, 
can be determined from the basic statistic, using the average value and the 
variance of the parameters of the battery output and so on which a 
designer cannot choose voluntarily. It was also confirmed that when at 
least either of the specified output of battery and the specified output of 
converter is determined taking the standard deviation into consideration, 
the system loss can be reduced. It was found to be particularly preferable 
that the specified output of battery is set to be in the range from not less 
than 0.7 time to not more than 2 times of standard deviation and the 
specified output of converter is set to be in the range from not less than 1 
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time to not more than 4 times of the standard deviation. This is because 
when the specified output of battery and the specified output of converter 
are in these ranges, respectively, the system loss is reduced further, 
leading to further improved system efficiency. 

In the system of this example, if variation is smoothed to be about 1.5 
times to about 2 times as in the example by e.g. increasing the number of 
windmills and time windows to be smoothed, then it is estimated that the 
converter loss will be reduced (5 — 2.5)/5 = about 50% (Cf. FIG. 4) and the 
battery loss will be reduced (15 — 8)/15 = about 47 % , thus achieving 
reduction of the system loss. In other words, in the system of this example, 
the converter loss+battery loss=5 + 15=about 20kW loss at 44.4kW which 
is the average of the absolute value of the battery output is reduced to 2.5 
+ 8 = about 10.5kW in total. Thus, the loss rate of the system can be 
reduced to approximately half, from 20/44.4=about 45% to 10.5/44.4=about 
24%. 

(Test sample) 

The system loss was examined, varying the specified output of the 
redox flow battery and the specified output of the DC/AC converter used in 
the example mentioned above. In the test, the power generation by wind 
was used for all samples, with its maximum output set at 400kW and its 
specified output of power generation set at 285kW. The specified outputs 
of the redox flow battery and converter used in the test, and the loss of the 
battery system used in the test are shown in TABLE 1. 
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TABLE 1 



Sample 


Specified 


Specified 


Loss of 


Smoothing of power 


No. 


output of 


output of 


battery 


generation by wind 




battery 


converter 


system 






(kW) 


(kW) 


(kW) 




1 


170 


275 


20 


Good 


2 


85 


147 


10 


Good 


3 


60 


200 


10 


Good 



In the samples No. 2 and No. 3, the specified output of the redox flow 
battery and the specified output of the DC/AC converter were determined 
5 from the standard deviation of the battery output determined in the 

example mentioned above. In the sample No. 2, the specified output of the 
battery was increased by 1.6 times of the standard deviation of 54.1kW and 
the specified output of the converter was increased by 2.7 times of the 
standard deviation of the same. In the sample No. 3, the specified output 

10 of the battery was increased by 1.1 times of the standard deviation of 

54.1kW and the specified output of the converter was increased by 3.7 
times of the standard deviation of the same. On the other hand, in the 
sample No. 1, the specified output of the battery and that of the converter 
were determined without any particular consideration of the standard 

15 deviation of the battery output. It was found from the test results that the 

samples No. 2 and No. 3 reduced the battery system loss to approximately 
half of the sample No. 1. 
Capabilities of Exploitation in Industry 

As explained above, the method of designing a redox flow battery 

20 system of the present invention can provide the advantageous result that 

the characteristic function to evaluate the system can be optimized by 
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determining an optimum value of the internal parameter based on the 
average value of the external parameter and the standard deviation. 

Particularly by determining the specified output of the battery and 
that of the DC/AC converter based on the average value of the battery 
output and the standard derivation, the system loss can be reduced. This 
enables the battery system loss to be reduced without annexing the lead 
storage battery to the system, as the related art does. This can provide a 
reduced scale of the hardware of the system and thus a more economic 
system, as compared with the related art. Further, since no lead storage 
battery is required, production costs can be reduced. 
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